Coil-helix transition in poly(L-glutamic acid) : Evidence for a 3-state non-cooperative 

process 
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A careful analysis of measurements of circular dichroism of poly(L-glutamic acid) (PGA) shows 
that the data can be very accurately described by introducing a third state for the PGA con- 
figuration, in addition to the helix and coil ones, and considering a simple equilibrium between 
these three states, without cooperativity. The third state is more conspicuous when high molecular 
weight polyethyleneglycol (PEG) is added. Excluded volume effects shown by differences in presence 
of short and long PEG chains indicate a direct interaction of PEG and PGA rather than an osmotic 
effect. 



PACS numbers: 87.15.Cc, 82.35.Jk 
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I. INTRODUCTION 

Proteins are important building blocks of life. They 
have multiple properties due to the interplay of their four 
structural levels : primary (sequence of aminoacids), sec- 
ondary ( alpha helices, beta sheets, etc.) tertiary ( ar- 
rangements of secondary structures in space) and quater- 
nary ( self assembly of supramolecular structures ) . The 
primary structure is fixed by covalent bonds and is there- 
fore very stable. The other levels of order are governed 
by much weaker interactions (hydrogen bonds, dipolar 
interactions, etc.) and are therefore sensitive to changes 
in the environment such as temperature, pH, etc. Under- 
standing the formation of these structures is an impor- 
tant step in the understanding of their overall properties. 
Moreover it is strongly suspected that misfoldings are re- 
sponsible for serious diseases [lj .We focus here on the for- 
mation of alpha helices from random coils in a much sim- 
plified system : the poly(L-glutamic acid) (PGA). The 
helical parts are optically active and their population can 
be easily and accurately measured by circular dichroism. 

This transition is usually analysed in terms of a theory 
due to Zimm and Bragg 0, [1] ( or variants thereof 0- 
6]). In this theory an elementary segment of the chain 
can exist in two states: either as as disordered chain or 
as a loop. A free energy difference between these states 
comprises energy and entropy terms : 



AF = AE — TAS 



(1) 



where AE and AS can be considered as independent 
of the temperature in a first approach. (A more complete 
description Q introduced a linear temperature variation 
in AE and a logarithmic one in AS*. The non linearity 
in the log function in the range from 280K to 370K is 
too small to be seen so that the equation applies with 
a different meaning of E and S). We choose the energy 
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FIG. 1: Measured dichroism (in degrees. cm 2 /mole) . From 
top to bottom : 0%, 5%, 7.5%, 10%, 20% and 30% PEG. Red 
lines are fits using the 3state non-cooperative model. Horizon- 
tal lines show the range of possible low-temperature asymp- 
tots (see text). 



unit such that fc# = 1. In the absence of other inter- 
actions, the relative ratio of states would be given by a 
Boltzmann factor exp(— AF/T). In addition, Zimm and 
Bragg introduced an interaction between neighbouring 
elements which facilitates the formation of a new loop 
beside an existing one and therefore the transition of the 
whole chain. The problem is therefore similar to the ID 
Ising model and can be solved. 

Many sets of data have been analysed using this model 
d, l3-fl0l| and significant cooperation parameters deter- 
mined even though a glance at the figures shows a rather 
poor fit to the data. A new set of measurements us- 
ing solutions with added polyethyleneglycol (PEG) and 
at pH=3.75 have been performed using state-of-the-art 
equipment [ll| and will be the basis of our analysis. The 
results are plotted in figure 1. 

Without PEG, the transition occurs at low tempera- 



III. 



THE 3- STATE MODEL 



10 1 F 



10" 







- 



10' 



10" 




2.8 



3.2 
1/T 



3.4 



3.6 

•10 -3 



FIG. 2: Same data as in fig 1 expressed as the ratio of optically 
active peptide to optically inactive ones vs the inverse tem- 
perature. From bottom to top: 0%, 5%, 7.5%, 10%, 20% and 
30% PEG. Red lines are fits using the 3state non-cooperative 
model. 



ture and only a part of it can be studied. The addition 
of PEG shifts the transition to higher temperatures but 
modifies the system. The figure suggests however that 
these modifications do not change qualitatively the basic 
mechanism. 



II. CRITICAL PLOT OF THE DATA 

A more clear picture can be obtained by plotting the 
data in a different way : Eql can be rewritten as : 



log[x/(l - a:)] = AE/T - AS 



(2) 



where x is the fraction of helices. In the absence of 
cooperativity a plot of log(x/(l — x)) vs 1/T would be 
a straight line and strong cooperativity would lead to a 
sigmoidal shape. 

Computing the helical fraction is, however, not 
straightforward because it requires the knowledge of the 
asymptotic values reached by the curves of fig 1 at high 
and low temperature, corresponding to fractions of he- 
lices of 0% and 100% . These are the optical properties 
of given molecular conformations and should not vary ap- 
preciably with temperature. The high temperature limit 
can be safely assessed for the data with 0% and 5% of 
PEG . A value of 37.4° has been published [H[ for the 
helix maximum contribution (difference between the up- 
per and lower asymptotes). We allowed for a variation 
of a few degrees of each of these values to get the most 
consistent pattern. 

We notice first that all the curves of fig 2 have a simi- 
lar shape: two linear parts at low and high temperature 
matched by a smooth crossover. 



Obviously the data cannot be fitted by any of the 
shapes predicted by the theory of Zimm and Bragg. How- 
ever, the presence of two straight asymptotes hints to the 
occurence of two transition free energies implying three 
thermodynamic states. The fraction of helices, in the 
absence of cooperativity, is easily computed from Boltz- 
mann factors : 



or 



x exp(-AFi/T) 
(1-x) ~ 1 + exp(-AF 2 /T) 



exp-(A£Ji/T- ASi) 



(1-x) l + exp-(AE 2 /T-AS 2 ) 



(3) 



(4) 



taking the coil state as the reference one.( The curvature 
of the lines implies that only one state is optically active.) 

Each curve can be fit very accurately using 6 inde- 
pendent parameters. The values found for the limiting 
values and the energies are reasonable. However, many 
of these values cannot be considered as precisely deter- 
mined. The high helicity side of the PEG poor samples 
and the low helicity side of the PEG rich samples are 
not sufficiently present in the experimentally accessible 
range. Moreover a slight change in a limit value leads 
to a change in the corresponding energy without altering 
significantly the quality of the fit. Among the many pos- 
sible ways of constraining the fit, we choosed to explore a 
possibility suggested by the figures namely constant en- 
ergies ( slopes) for all concentrations. The trial energies 
are taken from the best fit of the 'best balanced' curve at 
10% PEG. They are 10400/f and 7780^ or 21kCal.M~ 1 
and 15.5kCal.M~ 1 The fit of the experimental data with 
this formula can be considered as quite satisfactory for 
all the samples including PEG as can be seen in fig 1 and 
fig 2. The data without PEG however did exhibit a sys- 
tematic deviation and a different set of values of energies 
had to be found . The values are 13500^ and 8870if or 
27kCal.M~ 1 and 17.8fcCaZ.Af- 1 . This fit is the one plot- 
ted in fig 1 and fig 2. This 'constant energy' constraint 
is only a heuristic and fits with constant entropies and 
adjustable energies are as good. 



IV. EFFECT OF PEG CHAIN LENGTH 

The results of a similar study have been published a 
few years agoJToj- The main difference is that its au- 
thors used short PEG chains while Koutsioubas et al. 
ll| use long PEG chains . These data (picked from the 
published figure) can be plotted using the same process 
(figure 3). A look at figure 3 reveals similarities and dif- 
ferences. First the pattern of two transitions matched by 
a crossover is kept, but the third state is much less con- 
spicuous. We shall therefore not try to perform any quan- 
titative analysis for this state. The energy difference be- 
tween the coil and helix states (slope of the data at small 
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FIG. 3: Ratio of optically active peptide to optically inac- 
tive ones from Stanley and Strey[l(|. From bottom to top : 
10%, 15%, 20%, 25%, 30% and 35 % PEG ( of low molecular 
weight) 



1/T) is equal to that observed for long chains. A striking 
feature is that the addition of some short PEG has al- 
ways a roughly constant effect, while the shift of the data 
due to adding long PEG clearly levels-off. To illustrate 
this point we have plotted the temperature T\ — E\jS\ 
vs. the fraction of PEG for all the curves (figure 5) . The 
data cannot be compared one to one because of the dif- 
ferent pH but the trend of the curves seems significant. 
The idea is that a peptide can interact with more and 
more small molecules when the concentration of these 
molecules is increased, while the excluded volume effects 
prevent this for macromolecules. This shift has there- 
fore probably an entropic origin. A mechanism related 
to osmotic stress, as suggested by Stanley and Strey 
should depend on the concentration of PEG but not on 
the length of the chains. 



ones (32 or more residues). This clearly favors a third 
state induced by a self interaction of the chains. The 
characterisation of the third state by neutron or X-ray 
scattering can be envisionned but the simultaneous pres- 
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FIG. 4: Temperature Tl at which the ratio helix /coil = 1 
vs PEG concentration. Blue squares : long PEG chains , red 
triangles : short PEG chains 



ence of several conformations would make the unraveling 
of the spectra difficult. A X-ray study by Muroga et 
al [13| concluded that the observed data did not fit with 
what was expected for a coil when a coil was expected. 
Indeed the observed spectra decrease more slowly than 
expected at large wavevector indicating a more compact 
state of the scattering objects. A study of the dynamics 
of the helix-coil transition for PGA also concluded that 
this process was not a two-state process, [lij 



VI. CONCLUSION 



V. FINAL REMARKS 

The presence of a third state betweem coil and helix 
should not be too surprising. Indeed the amino acids 
involved are very prone to interacting with each other 
leading to different structures. A hairpin structure, pre- 
cursor of a beta sheet , could be a likely hypothesis. A 
similar study with different lengths of another peptide 
molecules (AEAAKA) has been published [1, m- Our 
analysis of the data from the published figures shows 
that the third state is not visible for the shortest chains 
(up to 26 residues) but seems to appear for the longer 



The first conclusion of this study is that the sharpen- 
ing of the transition due to first neighbour interaction as 
predicted by Zimm and Bragg Q is nowhere seen. In- 
stead the consideration of three states in mere chemical 
equilibrium describes accurately the data. The nature 
of the third state, which may have a significant impor- 
tance in health related problems remains to be investi- 
gated as well as the exact interaction of PGA with PEG. 
The evidence of excluded volume effects indicates a di- 
rect interaction. Circular dichroism measurements for 
PEG concentrations between 0% and 5% could be useful 
for this purpose. 
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